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Abstract 

We employ N-body/3-D gas dynamic simulations of the formation of galaxy clusters 
to determine whether cluster X-ray morphologies can be used as cosmological constraints. 
Confirming the analytic expectations of Richstone, Loeb, & Turner, we demonstrate that 
cluster evolution is sensitive to the cosmological model in which the clusters form. We 
further show that evolutionary differences are echoed in the gross morphological features of 
the cluster X-ray emission. 

We examine current -epoch X-ray images of models originating from the same initial 
density fields evolved in three different cosmologies: (i) an unbiased, low density universe 
with Oo = 0.2; (ii) an unbiased universe dominated by vacuum energy with Oo = 0.2 and 
Ao = 0.8 and (iii) a biased Einstein-deSitter model (0 = f , erg = 0.59). Using measures of 
X-ray morphology such as the axial ratio and centroid shifting, we demonstrate that clusters 
evolved in the two low Oo models are much more regular, spherically symmetric, and centrally 
condensed than clusters evolved in the Einstein-deSitter model. This morphology-cosmology 
connection, along with the availability of a large body of cluster X-ray observations, makes 
cluster X-ray morphology both a powerful and a practical cosmological discriminant. 

Subject Headings: cosmology: theory — galaxies: clustering — intergalactic medium 
X-rays: galaxies — hydrodynamics — methods: numerical 

1. Introduction 

Over the past fifteen years, observational studies have demonstrated that ~40% of clus- 
ters exhibit evidence for multiple components (Geller & Beers f982. Dressier & Shectman 
f988. West & Bothun f990, Jones & Eorman f990). Extensive theoretical studies have been 
carried out to explore the implications of the dynamical youth of some clusters and the 
details of cluster evolution (e.g., Cavaliere et al. f986, Evrard f990, Katz & White f993). 
In particular, it has long been suspected that a connection exists between the evolution of 
clusters and the cosmological density parameter (Gunn & Gott 1972; Geller 1984); recent 
work by Richstone, Loeb & Turner (1992) has explored this relationship in detail. By mod- 
elling cluster growth as the collapse of Gaussian distributed, spherical density perturbations, 
they have shown that clusters form much earlier in low density (Oo ~ 0.2) models than in 
= 1 models. This approach to cluster evolution necessarily omits a realistic treatment of 
the physical details of cluster formation, leaving questions regarding relaxation timescales, 
the effects of cluster growth on morphologies, and the effects of the surrounding large scale 
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structure. We attempt to answer these questions with N-body/3-D gas dynamic simulations 
of galaxy clusters within their cosmic environment. 

In this Letter^ we demonstrate a morphology-cosmology connection for X-ray clusters 
and discuss its use as a cosmological diagnostic. The basic idea is that clusters in low 
cosmologies are older and, hence, appear more relaxed and symmetric than their counterparts 
in an = 1 universe. Richstone et al. calculate the difference in median cluster ages between 

00 = 0.2 and = f to be ~ 0.3 Hq""^ . Because the crossing time for a 10 keV gas in the central 

1 Mpc of a cluster is 0.6 billion years, this age difference corresponds to many sound crossing 
times. The observable effects of this age difference should be most apparent in the X-ray 
emission from the intracluster medium (ICM), because the isotropic pressure tensor of the 
gas drives it to lie on equipotential surfaces on a few sound crossing times, unless interupted 
by a merger event. To quantify this effect, we perform a suite of numerical cluster simulations 
(§2) using three viable cosmological backgrounds: (i) an unbiased, open universe Oo = 0.2; 
(ii) an unbiased universe dominated by vacuum energy Oo = 0.2, Xg = 0.8 and (iii) a biased, 
Einstein-deSitter model = f, erg = 0.59. We analyse the cluster formation histories and 
X-ray morphologies in §3, and demonstrate that there are distinct differences in the cluster 
populations grown in high versus low Oo cosmologies. We discuss caveats in §4 and provide 
a summary in §5. 

2. Numerical Models 

In current theories, the nature of large-scale structure is controlled by a surprisingly 
small number of adjustable parameters. The present density parameter Oo = SttG/Jo/SHq 
and dimensionless cosmological constant Ao = Ac^/3IIq set the rate of change of scale factor 
with time, and thereby determine the growth rate of linear perturbations (Peebles 1980). 
The Hubble constant Hq sets the normalization and, typically, one invokes a non-zero cos- 
mological constant only to retain a flat metric if < f , so that Xg = I — fto- In addition, 
there are parameters required to describe the initial spectrum of linear perturbations. Ex- 
cept for models which invoke topological defects to generate perturbations ( e.g.^ strings 
or textures), the common assumption is that of a Gaussian random density field described 
by a power spectrum P{k) = with 6^ the Eourier transform of the linear density per- 

turbations. Eor reasonably smooth spectra, the spectrum in a finite range of wavelength 
centered on a scale of interest can be approximated as a power law with effective spectral 
index n^ff= din P{k)/ din k. Eor rich clusters of galaxies, the scale of interest is close to 
the fiducial normalization scale of 8 Mpc {h= Hq/IOO km s~"^ Mpc~"^). Besides n^ff, 
the other parameter required to describe the initial density field is erg, the present, linear 
extrapolated variance < (Sp/p)^ >i/2 within spheres of radius 8 Mpc. 

We generate a set of eight initial density fields from a standard cold dark matter (CDM) 
spectrum using a 32^ grid with Bertschinger's (1987) constrained technique. The constraint 
guarantees cluster formation in the simulation volume by requiring that the smoothed density 
field (smoothed with a Gaussian filter with Rf = 0.2L) evaluated at the box center be 2.5 — 5 
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times the rms fluctuation on the flltered scale. We produce two density flelds for each separate 
simulation volume; the simulation volume is a periodic cube of comoving length iv = 30, 40, 
50 or 60 Mpc on a side (unless otherwise stated, all quoted lengths assume h = 0.5). Each 
initial realization is evolved in the three cosmological backgrounds listed above, making a 
total of 24 runs. 

Density flelds are sampled from an unbiased CDM spectrum with ilh = 0.5. For the 

1/2 

= 1 models, we use the reduced normalization erg = 0.59 (erg = (i^Spj p)'^^ on the Sh~^ 
Mpc scale), equivalent to a bias parameter h = 1.7. This is the normalization inferred 
from the observed adundance of rich clusters (White, Efstathiou & Frenk 1993, Zabludoff 
& Geller 1993); an unbiased model favored by the COBE measurements produces too many 
rich clusters. We begin the = 1 simulations at an initial redshift = 9 and reduce the 
initial density flelds by a factor of 6(1 + Zi) = 17. Eor each run, two sets of 32^ particles are 
displaced from a regular grid using the Zeldovich approximation applied to the generated 
density held. One particle species represents the baryons and the other the dark matter. 
The particle distributions are evolved in time using the combined N-body/3-D gas dynamic 
scheme P3MSPH (Evrard 1988). Gravity, PdV work and shock heating are incorporated 
for the baryons. Radiative cooling is ignored. The minimum effective spatial resolution of 
the simulations is ~0.005iv, which varies from 150 to 300 kpc for the box lengths employed. 
We approximate power on scales larger than the simulation volume by allowing for a DC 
component (typical initial amplitude ~4%) in the density fluctuations; the volumes are then 
propagated in cosmologies with slightly higher background density. 

In order to make direct comparisons of individual clusters evolved within different cosmo- 
logical scenarios, we evolve each of the initial density flelds in all three cosmological models. 
The initial redshifts for the low Oo runs are determined by requiring that the overall linear 
growth factor be the same as that in the Einstein-deSitter model. Requiring erg = 1 (no 
biasing) at the present epoch in both low Oo models leads to an overall growth factor of 17 
and starting redshifts of z^ = 47.5 for the Oo = 0.2 model and z^ = 23.0 for the Oo = 0.2, 
Ao = 0.8 model. Reproducing the observed rich cluster abundance implies a normalization 
(Tg > 1 (White, Efstathiou & Erenk 1993), but because we are not aware of any compelling 
mechanism for making galaxies less clustered than the mass, we have adopted a conservative 
"what you see is what you get" approach in the low density models. We assume a baryon 
fraction Oj = 0.1 for a// the models, allowing for some amount of coUisionless matter (either 
dark matter or galaxies) in each cosmology. 

Our approach has the advantage that the merger history of a particular initial density 
held is similar for all cosmological models. The merger rate depends on and A, but 
the proto-clusters or cluster 'building blocks' are (to zeroeth order) independent of the 
background cosmology. The price paid for this advantage is that the low Oo initial conditions 
are no longer consistent with a CDM spectrum. This is because features in the CDM power 
spectrum scale as the product 0/i, which is flxed to the value 0.5. The appropriate value 
for the low Oo runs is ilh ~ 0.2. Using ilh = 0.5 causes the local spectral index n^ff 
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Figure 1: Contour maps of cluster X-ray emission: Tlie first tliree columns show simulated clusters evolved in the 
following models: (i) Qo = 0.2; (ii) Qo = 0.2, Ao = 0.8; and (iii) Q = f, erg = 0.59. The fourth column consists of 
Einstein Imaging Proportional Counter (IPC) observations of X-ray clusters. Each row in the first three columns 
corresponds to the same initial density field. The contours in all sixteen maps are spaced by factors of f .8 in surface 
brightness from 3.5 x fO~* cts s~^ arcmin"^, and the spatial scale for every map is identical (the distance between 
tic-marks is f88 kpc). The IPC observations have been corrected for instrument response and Gaussian smoothed; 
from top to bottom, the clusters are: (i) A496, (ii) A399, (iii) A2256, (iv) and A40f. 

{n^ff = din P{k)l din k) to be slightly larger than it should be, roughly —1 rather than 
— 1.3. This local spectral index controls the shape of the density distribution in collapsed 
halos, but the difference in halo profiles for the above spectral indices is very small (Crone, 
Evrard & Richstone 1993). We note that this range of spectral index is consistent with the 
power spectrum slope inferred from large-scale structure in the Cf A redshift survey (Vogeley 
et al. 1992) and the 1.2 Jy and QDOT 5 surveys (Fisher et al. 1993; Feldman, Kaiser 
& Peacock). 

3. Comparison of X— ray Morphology 

From the gas particle data for each simulation, we calculate X-ray maps (projected 
along three orthogonal axes) in the Einstein 0.5 — 3 keV energy band following the procedure 
described in Fvrard (1990). No Poisson noise or observational artifacts are added. Figure 
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1 contains contour maps of the X-ray emission from half of the simulated clusters. The 
first three images in each row represent the same initial density field evolved to the present 
epoch in the three different cosmological models. Clusters evolved in the low Oo universes are 
generally more spherical and centrally condensed than those evolved in the Einstein-deSitter 
model. This difference results from the accelerated evolutionary histories experienced by 
clusters in the low density universes. 
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Figure 2: Merger hierarchy of the same initial density field (50a) evolved in the three different cosmological models. 
This plot follows the present-epoch cluster mass backward in time, identifying the protoclusters which eventually 
merge to form the final cluster. The vertical lines correspond to mergers of subclumps with the main cluster. The 
third row of Figure 1 shows contour maps of clusters formed from this initial density field. 

We demonstrate this effect explicitly for the clusters in Figure 2 which were evolved 
in a 50 Mpc box. In this merger history diagram, we follow the mass associated with the 
final cluster backward in time, identifying the protoclusters which merge along the way to 
form the final cluster (the vertical lines correspond to mergers of subclumps with the main 
cluster; the algorithm is described fully in Evrard, Summers & Davis 1993). The cluster 
evolved in the Oo = 0.2 model grows rapidly at early times (50% of its final mass is in place 
by tjto = 0.3) and relatively slowly at late times. In the Einstein-deSitter model, growth 
begins more slowly and continues into the present epoch; the cluster reaches 50% of its final 
mass at tjto = 0.7. Both the X-ray morphologies (Eigure I, row 2) and merger histories 
(Eigure 2) indicate that the Oo = 0.2, Xq = 0.8 model is an intermediate case which is more 
similar to the Oo = 0.2 model; the cluster contains 50% of its final mass at t/to = 0.4 — 0.5. 

Also shown in Eigure I are contour maps of Einstein X-ray images of four Abell clusters: 
A496, A399, A2256 and A40I. Comparing these to simulated clusters indicates that the 
observations better match the clusters evolved in the Einstein-deSitter model. To compare 
clusters in these three models more quantitatively, we characterize the morphology in terms 
of two parameters: the axial ratio and the centroid shifting (Eigure 3). Applying a method 
described in Mohr, Eabricant, & Geller (1993b), we allow for variations in centroid and axial 
ratio as a function of radius in the cluster by examining the images within annuli of constant 
width (I') and varying radius (from I' to 16'). We determine the correct location of an 
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annulus by minimizing the difference between the annuius center and the centroid of the 
photon distribution within the annuius. This approach avoids biasing the resuits by an a 
priori choice of image center. For each annuius, we then determine the image centroid and 
axial ratio. The axial ratio plotted in Figure 3 is the emission weighted mean axial ratio 
for the entire image, and A, the centroid shift, is the emission weighted rms variation in 
centroid for the entire image. These measurements indicate that the two low Oo models 
produce clusters which are more regular and spherically symmetric than those evolved in the 
Finstein-deSitter model. 

Figure 3: The distribution of ax- 
ial ratio (left column) and centroid 
shift (right column) for all models 
(§3): Each row corresponds to a 
different cosmological model. A, 
the centroid shift, is the emission 
weighted rms variation in centroid 
given in arc minutes (f ' = 70 kpc); 
the axial ratio presented here is 
the emission weighted mean axial 
ratio. An inner region with a 6' ra- 
dius is excluded in calculating the 
axial ratio and centroid shift of the 
fio = 0.2 models. 
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We are expanding on the work presented here by completing a systematic comparison of 
this battery of simulations with a representative ensemble of clusters observed by Einstein . 
By introducing the appropriate instrument reponse, accounting for X-ray point sources, and 
making the same measurements on simulated and observed clusters, we find that the clusters 
evolved in the Finstein-deSitter model most closely resemble observed clusters (Mohr et al. 
1993a). 



4. Discussion 

Simulations, by nature, provide very specific information. One immediate and important 
question is how generic are the above results? What are the key controlling elements and 
how might they affect the results? We explore these issues here. 

The rounder X-ray isophotes in the low Oo models result primarily from the scarcity 
of recent mergers; this allows the gas to relax onto isopotential surfaces. In addition, the 
isopotential contours are rounder in the low Oo clusters because the gas represents half of 
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the mass. In contrast, the = 1 clusters have 90% ol their mass in the form ol dark matter 
supported by an anisotropic velocity-dispersion tensor. 

There are several potential concerns arising from the finite simulation volume. One is 
that the low Oo clusters may have had their accretion shut ofi' artificially early. There are two 
lines of evidence discounting this idea. First, analytic models predict the early formation 
of clusters in low Oo universes (Richstone et al. 1992; Lacey & Cole 1993; Kauffmann 
& White 1993) and, second, the merger histories from these simulations agree well with 
merger histories calculated from an ensemble of clusters obtained from N-body simulations 
of larger volumes (Crone, Richstone, & Evrard 1993). This result is insensitive to the shape 
or normalization of the fiuctuation spectrum (Lacey & Cole 1993) for reasonable definitions 
of the cluster population. 

Another concern is that tidal effects generated on scales comparable to the simulation 
volume could distort the cluster gas, altering its X-ray morphology. However, extremely 
strong (and unrealistic) torques are required to affect gas in highly non-linear. X-ray emitting 
regions. For example, consider a point mass perturber of mass Mp a distance R from a cluster 
of mass Mf-i- For the gravitational tides to have a fractional effect / with respect to the cluster 
self-gravity on a test particle a distance r from the cluster center would require 



For values r = l Mpc and R = L/4: = 7 .5 Mpc for our smallest box, a 10% effect on isophotes 
at 1 Mpc would require Mp > 40Afci, which is unlikely. 

Missing physics could also affect our results. It is difficult, however, to come up with 
mechanisms capable of strongly perturbing the low Oo gas morphology away from spheri- 
cal symmmetry at the present epoch. Clusters evolved in the low Oo models have central 
densities sufficiently high that radiative cooling times are much less than the Hubble time. 
If radiative cooling (but no star formation and associated feedback) were turned on in the 
models, a large puddle of cold gas would collect in the cluster center, leaving a rarefied 
envelope of hot gas. There is no reason to suspect that the morphology of the hot envelope 
would be any less spherical than in the present models. Inclusion of galactic winds could add 
energy and metals to the gas but, again, generating an aspherical morphology would require 
(i) a directed energy input comparable to the total energy in the X-ray emitting gas (~ 10^^ 
ergs) and (ii) that galaxy formation occur after cluster formation, so that the feedback occur 
very recently. These conditions are extremely unlikely. Simulations incorporating realistic 
feedback from galactic winds in an = 1 model show there is little effect on the gross 
morphological features of the X-ray emission at the present epoch (Metzler & Fvrard 1993). 

Finally, we have not explored a wide range of values of Oo and so we cannot explicitly 
state the sensitivity of cluster morphology to the exact value of Oo- As Oo increases from 
0.2 toward unity, the linear growth rates rapidly approach that for the = 1 model. In 
addition, adding a Xg component to fiatten a model serves to make the growth curve more 
similar to the = 1 model. Our results, therefore, are better interpreted as an argument 
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against Oo = 0.2 than an argument for = 1. We hope to explore the Oo dependency of 
cluster morphology more extensively in future work. 

5. Summary 

We demonstrate a morphology-cosmology connection for X-ray clusters: evolutionary 
differences controlled by the underlying cosmological background are reflected in gross mor- 
phological features of the X-ray emission. This confirms the analytic expectation of Rich- 
stone, Loeb & Turner (1992), whose analysis was based on a spherical model of gravitational 
collapse. Using a combined N-body/3D gas dynamic method, we evolve eight different ini- 
tial density fields to the present epoch in three different cosmological models: (i) Oo = 0.2, 
(Tg = 1; (ii) Oo = 0.2, Xg = 0.8, erg = 1 and (iii) = 1, erg = 0.59. We determine the axial 
ratio and centroid shifting of idealized 'observations' of these clusters and use these param- 
eters to characterize the morphological differences in these models. Clusters evolved in the 
two low Oo models appear far more regular, spherically symmetric, and centrally condensed 
than clusters evolved in the Einstein-deSitter model. We are currently completing a direct 
comparison of simulated and real clusters using an ensemble of Einstein observations; results 
indicate that, of the three models tested, clusters evolved in the Einstein-deSitter model best 
match the complex nature of real cluster X-ray morphologies (Mohr et al. 1993a). 
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